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ABSTRACT 

Efforts to benchmark; astrophysical observations with X-ray laboratory measurements 
have been stymied by observed and measured differences of up to a factor of two in 
the ratio '3s/3d' of Fexvii Unes at AA ~ 17 A and ^ 15A respectively. Using the 
electron distribution function (EDF) as a new physical parameter, we compute the 
Fexvii line ratios and account for these differences. Based on large-scale relativistic 
close coupling calculations using the Breit-Pauli R-matrix method, revealing the pre- 
cise effect of resonances in coUisional excitation, we employ collisional-radiative models 
using cross sections convolved with both the Gaussian and the Maxwellian EDF. Com- 
parison with astrophysical observations and laboratory measurements demonstrates 
that (a) the 3s/3d line ratio depends not only on the EDF but also on the electron 
temperature/energy of the source, (b) plasma conditions in experimental measure- 
ments and astrophysical observations may be quite different, and (c) departure from 
a Maxwellian should manifest itself in, and be used as a diagnostics of, particle distri- 
butions in plasmas. 

Key words: Gaseous Nebulae - Optical Spectra: H ii Regions - Line Ratios: Atomic 
Processes - Atomic Data 



1 INTRODUCTION 



Astrophysical X-ray sources cover a wide range of physi- 
cal conditions in objects as diverse as stellar coronae, ac- 
tive galactic nuclei. X-ray binaries, supernovae, and after- 
glows from gamma-ray bursts. In principle, spectral diagnos- 
tics should provide accurate information on the particular 
plasma conditions, particularly from high-resolution obser- 
vations p ossible with X-ray sa t ellites Chandra and XMM- 
Newton ijCaim ares etal. 2000l : iHuenenmoerder etal. 20011 : 
ISaba etal. 19991') . However, the diagnostics are contingent 
upon accurate theoretical models calibrated against pre- 
cise experimental measurements. Several spectral lines from 
neon- like Fexvii are prominent constituents in the soft 
X-ray region ~15-17 A, and have long been observed as 
potential diagnostics (Fig. 1). Recent electron beam ion 
trap (EBIT) experiments on Fexvii line intensity ratios 
at the Lawrence Liverm ore National Laboratory (LLNL) 
jBeiersdorfer etal. 20o3) . and at the National Institute for 
Standards and Technology (NIST) jLaming etal. 200Cl) . 
may establish these diagnostics, provided the experimental 
results, theoretical calculations, and astrophysical observa- 
tions can be interpreted consistently. But the diagnostic ca- 
pabilities of the Fexvii X-ray line ratios have so far been 
limited because the line-formation mechanism is not well 
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Figure 1. Schematic Fcxvil transitions for three 3s X-ray lines 
(3F, 3G, and 3H) and three 3d X-ray Unes (3C, 3D, and 3E). 



understood, and all spectral modeling codes carry large un- 
certainties. 

A schematic representation of the lines of interest in 
this Letter is given in Fig.0 Excitation of the closed-shell 
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2p^ ground configuration of neon- like Fe xvil to tlie 2p^ 3s 
and 2p® 3rf configurations corresponds to sets of 3 lines 
each at ~ llA and ~ 15A, respectively labeled as the 
'3s' and the '3d'. The interpretation of the 3s/3d ratios 
in Fexvil from current EBIT and tokamak experiments, 
and astrophysical observations, are perplexing. Basically, 
there are two long-standing problems in Fexvil spectra: 
the so called 3C/3D ratio and the 3s/3d ratio in Fexvil 
jLaming etal 2000l: ISeiersdorfer etal. 2002^ . The 3C/3D is 
the line intensity ratio of transitions labeled as 3C (A 15.015 

A) and 3D (A 15.262 A). The ratio 3s/3d=i|^±|^±||l 

(line 3H is sometimes labeled as M2), is the intensity ratio 
of three lines at AA 17.098, 17.053 and 16.778 A with one 3s 
electron in upper levels, over another set of three lines at AA 
15.015, 15.262 and 15.450 A with one 3d electron in upper 
levels. 

In the past theoretical predictions of the X-ray line 
intensity ratio 3C/3D were about 50% or more higher 
than astronomical observations or laboratory experiments. 
Yet, the EBIT measurements were in agreement with each 
other and with observations. Th is problem has been ad- 
dressed in our recent calculat ions dChen and Pradhan 20o3 : 
IChen etal. 20031 : IChen 2005^ in the low electron impact en- 
ergy region. In this low energy region, we carried out ab 
initio cal culations using the Bre it-Pauli R-matrix (BPRM) 
method ijBerrington etal. 1995^ including electron-impact 
excitation (EIE) resonances up to n = 4 levels (n is the 
principal quantum number) , which play a key role in the so- 
lution of the 3C/3D problem. The resonance enhancement 
is extrapolated from n = 4 to n = 5 for the higher electron- 
energy (temperature) region (this extrapolation needs to 
be verified by further calculations). Similarly, the 3s/3d 
problem pertains to nearly a factor of 2 difference in the 
intensit y ratio between the two sets of EBIT m easure- 
ments jBeiersdorfer etal. 20021 : iLamine et al. 200(11 . More- 
over, there is very large scatter in the 3s/3d ratios from all 
previous theoretical calculations and observations. However, 
this case is more complicated than the 3C/3D problem for 
two reasons. First, at low electron densities cascade effects 
have important or dominant contributions to the formation 
of some or all 3s X-ray lines; while the 3C and 3D lines are 
formed predominantly by EIE. Second, the 3s/3d ratio is 
very sensitive to plasma conditions under which these X- 
ray lines are formed. In contrast, the 3C/3D ratio is not 
very source sensitive i.e. it is determined mainly by atomic 
micro-processes. It is due to these two reasons that the mea- 
sured and observed 3s/3d ratios for given temperature and 
density may differ by up to a factor of two — a difference 
that is a major unresolved problem in X-ray astrophysics. 

One of main goals in this paper is also to introduce 
electron distribution functions (EDF) (e. g. Maxwellian and 
Gaussian EDF) as a new physical parameter for the study of 
X-ray formation mechanism. As evident from the results in 
this Letter, this new approach is essential to solve the long- 
standing 3s/3d problem in Fe XVII. The main point is that 
plasmas in EBITs and astrophysical sources may have dif- 
ferent beam/particle distributions. In some laboratory plas- 
mas e. g. EBIT experiments, the electron beam width may 
be Gaussian or non-Gaussian (non-Maxwellian) . In astro- 
physical objects on the other hand the EDF is normally a 
Maxwellian, but there are important cases (such as mag- 



netically driven plasmas) where the source may display a 
bi-Maxwellian character, or a non-Maxwellian component 
of electron fluxes with high-energy tails. The precise con- 
volution of cross sections with velocity distributions must 
therefore be generally source-specific. 



2 THEORY AND COMPUTATIONS 

The extensive calculations of the Fe xvil collision strengths 
using the BPRM method and a 89-level close coupling 
expansion up to the n = 4 levels is described in ear- 
lier works dChen and Pradhan 20o3lChen etal. 2003h . Res- 
onance structures in the collision strengths are delineated; 
the positions, strengths, and distribution of resonances have 
been illustrated in detail for many of the transitions of inter- 
est in this work. Considerable differences are found with re- 
spect to previous calculations that neglect resonances, such 
as several distorted wave results available in literature. The 
energy dependence of the collision strengths due to reso- 
nances introduces a temperature dependence in the aver- 
aged collision strengths with temperature dependent EDF's 
such as a Maxwellian. Moreover, the many Rydberg series 
of resonances, converging on to most of the 89 levels in the 
wavefunction expansion, are not uniformly distributed in en- 
ergy. Therefore a Gaussian average yields significantly differ- 
ent values, depending on the assumed energy widths, than 
a Maxwellian average. 

In the present work, in order to address the issue of EDF 
influence on line ratios in general, we recalculate the aver- 
aged collision strengths, or rate coefficients, and reconstruct 
the coUisional-radiative model. We carry out separate convo- 
lutions of our detailed close-coupling cross sections with spe- 
cific beam distributions e. g. a Gaussian average (GA) and 
a Maxwellian average (MA) , before the construction of the 
coUisional-radiative model. The Maxwellian averaged rate 
coefficients are computed for all transitions as 



8.63 X lO^'^cm^/s • exp{-Eij/kTc) 



(1) 



where gi is the statistical weight of the initial level and the 
quantity Tij is the Maxwellian averaged collision strength: 



T,,(T) 



n,j(S)exp(-S/fcTo) d{E/kT,). 



(2) 



In addition, calculation for all rate coefficients are re- 
peated with Gaussian EDF's obtained using a specified en- 
ergy spread AE to simulate different beam widths in EBIT 
experiments. The coUisional-radiative model for the Fexvil 
lines includes 3,916 transitions associated with the atomic 
model, which also employs new allowed and forbidden tran- 
sition probabilities for Fe xvil recently computed by Nahar 
etal. (2003). 



3 RESULTS AND DISCUSSION 

The main results of the calculations and modeling of the 
3s/3d line ratio are summarised in Fig.|21 The black arrow 
marks Tm ~ 4 x 10® K, the value of Fexvil maximum abun- 
dance in coronal equilibrium. The red arrows mark different 
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highest target thresholds that may be included in CC calcu- 
lations. Although this study is based on one of the largest 
and most computer intensive BPRM close coupling calcula- 
tions, we could only include target levels up to n = 4. Com- 
parison with the two experiments shows that NIST EBIT 
data mimic our calculation assuming a Maxwellian averated 
EDF (MA), and LLNL data are in basic agreement with 
the GA results. However, there is some discrepancy with the 
LLNL results at the higher temperature end. The target lev- 
els with 71 ^ 5 are currently neglected in the 89CC BPRM 
calculations. This limitation could result in some CEiscad e 
contribution to the 3s/3d X-ray line ratios ijLiedahl 199fll . 
as well as some missing resonance structures. We have used 
an EBIT electron beam distribution with FWHM (full width 
at half maximum) of 30 eV in our GA modeling. 

Fig. 121 shows two sets of 3s/3d intensity ratios as func- 
tions of electron temperature or electron beam energy. The 
first set of calculations using MA collision strengths in the 
89-level CR model appear as the black curve, calculated at 
two electron densities 10^"^ cm"'' (black solid line) and 10^ 
cm~^ (black dotted line); dependence upon electron density 
is negligible. The second set of data are shown by the blue 
solid curve with 3s/3d ratios using GA collision strengths 
with a FWHM of 30 eV (blue dotted curve is obtained with 
a 50 eV convolution which gives essentially the same re- 
sults); the further calculation with a 5 eV convolution is 
also shown by the red solid curve. The oscillations in the 
3s/3d values in the GA model are clearly seen. They are 
due to the distribution, strength, and the density of reso- 
nances in the collision strengths in different energy regions 
as sampled by the convolution function centered around a 
30 eV or 5 eV FWHM Gaussian. This oscillation behavior 
is crucial in explaining accurately the energy (temperature) 
dependence of 3s/3d ratios in the EBIT experiments. All 
previous calculations either fail to get the correct 3s/3d line 
ratios, or fail to obtain their right temperature dependence. 
The observed 3s/3d ratios from the solar corona are marked 
with red circles. The open and filled blue circles with error 
bars mark in particular the coronal spectra of the solar-type 
star Capella ^ observed from Chandra and XMM-N ewton 
respectively iCanizares etal. 2000l : IXudard etal. 200 ll) . The 
LLNL EBIT experimental values are shown as green filled 
circles with error bars, and the NIST EBIT experimen- 
tal values as tw o black open circl es (without polariza- 
tion corrections jLaming etal. 2000l) }_ or two green op en 
circles (with polariz ation corrections l|Gillaspv etal. 20041) 1. 
N. B. NIST data in (iLaming et al. 2000l) are reported with- 
out polarization correct ions while all the LLNL data in 
jBeiersdorfer etal. 20o3) are reported with polarization cor- 
rections. Also shown for comparison are the very scattered 
3s/3d ratios from the long history of various observations of 
solar, stellar, and disk coronae (corresponding references are 
given in the caption of Fig. 2) . 

The 3s/3d problem arises mainly from the factor two 
difference in the two EBIT experiments at the lower en- 
ergy, and partly from the large scatter in other observa- 
tions. This problem is more serious than the 3C/3D prob- 
lem, as mentioned above, because the two EBIT group are 
in good agreement for the 3C/3D ratio. Also, the 3s/3d ra- 
tio is complicated by a complex line formation mechanism 
in the three 3s X-ray lines. Both, the cascade effects from 
high-lying levels and resonance effects, are important and 



need to be calculated accurately. First of all, we need to 
understand the difference in EBIT experimental values. We 
note that the 3s/3d values from EBIT measurements need 
to be interpreted with care in order to compare directly 
with astronomical observations because of different plasma 
conditions and/or different EDF. This point does not bear 
strongly on the 3C/3D problem, but for the 3s/3d prob- 
lem it is critical in order to derive a correct interp retation. 
For example, a recent work (iBeiersdorfer etal. 200411 reports 
similar Y&hies. for the line ratio 3s/3C (which differs slightly 
from the ratio 3s/3d in the denominator) from two differ- 
ent sources, the LLNL EBIT and t he tokamak Princeton 
Large Torus iBeiersdorfer etal. 2004[). However, this kind o f 
agreement as detailed in Fig. 7 of iBeiersdorfer et al. 2004^ 
is most likely fortuitous. This difference in magnitude for 
GA results, decreasing monotonically with temperature over 
extended ranges but attenuated by oscillations due to res- 
onances in specific energy regions, is clearly seen in Fig.|21 
and is the general feature established in this work. At lower 
energies the effect of resonances belonging to low-lying lev- 
els is more pronounced, and decreases with successively 
higher levels. This asymmetry in the distribution of reso- 
nance strengths with energy is responsible for the overall 
behavior of the line ratio in Fig.|21 

Fig-IHapparently shows that the NIST data at the lower 
energy agrees better with the primary MA calculations, al- 
though the EDF inside both the NIST and the LLNL EBITs 
is quite different from a Maxwellian. The dichotomy inher- 
ent in these results reveals the importance of precise deter- 
mination of experimental beam energies and widths. The 
NIST data point at the lower energy could possibly corre- 
spond to the minimum in the fundamental oscillations of the 
3s/3d line ratio, with a narrower EDF than the LLNL EBIT 
iCillaspv etal. 2004ll . We further illustrate this by calculat- 
ing a GA 3s/3d line ratio with a FWHM of 5 eV as shown 
in the oscillating red solid curve in Fig.|21 We find that in- 
deed the NIST data (the green open circle with polariza- 
tion correction) at the lower energy is located very close to 
a minimum in this red curve. The NIST data point (with 
polarization correction) at the higher energy is in agree- 
ment with our GA result. While at present the two sets 
of experimental data do not fully trace the theoretically 
predicted oscillations, the present results show that a much 
more detailed investigation of the fundamental nature of the 
measurements is required in order to reconcile the two sets 
of EBIT measurements, assuming that both are of claimed 
accuracy. Accurate measurements of beam widths are cru- 
cial, but to our knowledge have not been accurately deter- 
mined or reported. The EDF, in turn, bears on physical 
effects such as space cha rge potential, beam currents, polar- 
ization corrections etc. iCillaspv etal. 2004^ . Nevertheless, 
the theoretical sets of 3s/3d ratios clearly demonstrate that 
they are sensitive to the pl a sma conditions in the s ource 
(c.f. (IBeiersdorfer etal. 20021: IBeiersdorfer etal. 2004) 1. In 
order to interpret additional experimental information on 
these issues, the theoretical results will also need to be ex- 
tended to other EDF to elicit a more complete description 
of resulting variations in line intensities. 

The large-scale relativistic CC calculations for EIE 
of X-ray lines in Fexvil llChen and Pradhan 20021 : 
IChen et al. 2003^ . and the study of different electron dis- 
tribution functions in plasmas via CR models, appear 
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Figure 2. Theoretical X-ray line ratios 3s/3d as a function of electron temperature, or electron beam energy E = 
kT, are com pared with observatio ns from the corona of solar- type star C apella at ^ 5—6 MK from Chandra (open 
blue circle) iCanizares etal. 20001) and XMM-Newton (filled blue circle) jAudard etal. 200 J) , from solar coronae 
at Tm ~ 4 MK (red open circles) jHutcheon etal. 197fi). and from the LLNL and NIST EBIT experiments (filled 
and open green circles respectively') iBeierado^T etal. 2002HLaming etal. 200Cl) ^. Sec text for more discussions 
about the details of the 3s/3d problem. The scattered data for 3s /3d line ratios wit h labels 1-14 are for various 
types of observed values from solar, stellar and disk c oronae tMckenzie_efaLJ^80|; |Huenenmoeixieii_efaL_2001 ; 
Ness etal. 2003tlMewe etal. 200]l:lBrinkman etal. 2000l: iRaassen etal. 2002HPhilliDS etal. 1997HKahn etal. 2001 ; 
Xu etal. 2002^ . 



to have resolved long-standing problems. The relativistic 
atomic calculations of Fexvil and the line ratios among 
its diagnostic X-ray lines are in agreement vyith two in- 
dependent measure ments on EBIT (iBeiersdorfer etal 20o3: 
iLaming etal. 2000l) . It follows that the line ratios from the 
present 89CC CR model may be further applied to interpret 
X-ray spectra from a variety of X-ray sources observed from 
current space observatories, Chandra and XMM-Newton, 
and the upcoming next generation ones with greater spec- 
troscopic capabilities, such as Constellation-X. Towards this 
goal additional calculations are in progress to calculate level- 
specific and total recombination rates from from Fe XVIII 
to Fexvil , including radia tive and dielectronic recom bina- 
tion in a unified scheme (iNahar and Pradhan 20(3) . and 
photoionization rates from Fe XVI to Fe XVII. However, 
astrophysical environments and plasma conditions may be 
different from those in EBIT plasmas. It is our goal to ap- 
ply our 'calibrated' atomic data to interpret observed stellar 
and cosmic X-ray spectra of coronal or coUisionally ionized 
plasmas such as stellar coronae, and photoionized or hybrid 
plasmas in, for example, from active galactic nuclei (AGN), 
quasars, disk coronae etc. 



In summary, the new results for Fexvil structure, col- 
lision dynamics, and spectral diagnostics reveal not only the 
effect of resonances and cascades, but also present a chal- 
lenge for the study of source-specific effects of electron dis- 
tributions and other conditions prevalent in EBIT, tokamak, 
laser-produced, and astrophysical plasmas. The present GA 
and MA averaged CR models should offer possible solutions 
to obtain more precise information on the physical param- 
eters in the beam and the trap in EBIT experiments. We 
demonstrate that 3s/3d X-ray line ratios are not a univer- 
sal constant in an optically thin plasma excited by electron 
impact, but instead depends sensitively on both the charac- 
teristic electron energy (temperature) and the detailed dis- 
tribution of the electron energies about that characteristic 
value. Furthermore, the present results should be of general 
astrophysical importance for objects in a variety of equilib- 
rium or non-equilibrium conditions. A direct physical conse- 
quence of this work is that the departure from a Maxwellian 
should manifest itself in, and be used as a diagnostics of, 
particle distributions in plasmas. In photoionization equi- 
librium (as opposed to coronal) we also need level-specific 
recombination cross sections and rates in order to compute 
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emission line intensities (work in progress). Finally, the re- 
sults in this Letter should be appli cable to other ne on-like 
ions of importance in X-ray lasers ijRosen etal. 198^ . with 
population inversion among the upper n=3 levels due to EIE 
and/or electron- ion recombination. 



4 CONCLUSION 

The main conclusions of the present study are: 

• Line intensities and ratios of the 3s and 3d line- 
complexes at ~ 17 A and ~ 15A depend on electron distribu- 
tion functions. The 3s/3d ratio can be a potential diagnos- 
tics of temperature, particularly for T ^ Tm in astrophysical 
sources with Maxellian EDF's, and of energy, particularly 
at i5 > kTm in non-Maxwellian plasmas such as laboratory 
sources. 

• The Maxwellian and Gaussian simulations of Fe xvil 
X-ray line ratios provide an indication of diffrerences in the 
conditions in the interaction regions of Electron-Beam-Ion- 
Traps; however, more experiments are needed to fully un- 
derstand the physical factors involved. 

• X-ray line ratios are a potential diagnostic of non- 
maxwellian astrophysical and laboratory sourcs, such as ac- 
cretion discs. X-ray flares and bursts, and magnetically or 
inertially confined fusion plasmas. 
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